Electron tunneling through graphene-based double barriers driven by a periodic 

potential 
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Photon-assisted charge transport through a double barrier structure under a time periodic field 
in graphene is studied. Within the framework of the Floquet formalism and using the transfer 
matrix method, the transmission probabilities for the central band and sidebands are calculated. A 
critical phase difference between the harmonic potentials at the barriers, which cancels transmission 
through the inelastic sidebands due to quantum interference is found. This phenomenon could be of 
help to design graphene based filters and high-frequency radiation detectors. Quenching of resonant 
tunneling by the harmonic field applied to the barriers or the well is also discussed. 



Quantum transport in periodically driven mesoscopic 
systems is an important subject not only of academic 
value but also for device applications. The interest on 
time-dependent excitations by electromagnetic fields has 
been increasing and many interesting phenomena have 
been investigated [l|, 0. Early studies of Dayem and 
Martin Q provided the evidence of photon assisted tun- 
neling (PAT) in experiments on superconducting films 
under microwave fields, and subsequently Tien and Gor- 
don theoretically justified this observation After 
this, electron transport through various types of time- 
oscillating potential regions has been studied in semicon- 
ductor nanostructures [5H8| aiming to the fabrication of 
new devices. 

A few years ago, graphene, a two dimensional honeycomb 
crystal of carbon atoms, was fabricated by Novoselov et 
al. @, leading to increasing attention to both the fun- 
damental physics and potential applications of the new 
material. In graphene, low-energy quasiparticle exci- 
tations near the edges of the hexagonal Brillouin zone 
(Dirac points) obey a gapless linear dispersion law, and 
their motion can be described by a two-dimensional Dirac 
equation for masless particles. The presence of such 
Dirac-like quasiparticles leads to Klein tunneling and 
other unusual electronic properties ffjl-[l4j . For graphene- 
based PAT devices it is essential to consider transport 
of charge carriers in graphene through time-harmonic 
potentials. PAT through single barrier in monolayer 
graphene have been discussed in Ref. This letter 

is an extension for double barrier structures (DBS) tak- 
ing into account quantum interference between photon 
assisted processes. At normal incidence PAT is analyzed 
aiming to get a direct insight into the quantum inter- 
ference and analytical condition for the suppression of 
the transmission through the inelastic sidebands is ob- 
tained. We also find that quantum interference makes 
a significant contribution to the total transmission prob- 
ability at non-normal incidence and could be useful in 
graphene-based PAT devices. 

A single electron transmitting through a monolayer 
graphene-based double barrier driven by a harmonic po- 
tential is considered. This structure can be fabricated by 



FIG. 1: Schematic oscillating potential of the DBS in 
graphene. 



applying a local top gate voltage and a small ac signal to 
graphene. At low energy, and close to the Dirac points 
(K and K') electrons are described by a massless Dirac 
Hamiltonian: 



with 



H = Vfd ■ p + V(x, t) 



V(x, t) = I [V + V{x)cos (ut + 6>d (x) )] 



(1) 
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where the Fermi velocity Vf « 10 6 m/s, a — (a x ,<j y ) is 
the 2D vector of Pauli matrices and p — > —ihV is the 
2D momentum operator. The barriers of height Vq are 
oscillating with frequency u, phase difference S, and dif- 
ferent amplitude: V(x) = Vt if x < or V[x) — V3 for 
x > 0. This simple device is shown schematically in Fig. 
1. The barriers and well widths are L\ = (L — d)/2 and 
L-2 = d respectively, $(x) is the Heaviside function and 
I is unit matrix. The DBS is infinite and homogeneous 
along the y-direction, resulting in the conservation of the 
y-component of momentum. Incident electrons propa- 
gate from left and pass five regions denoted by 0, 1, 2, 3, 
4. The solution of Dirac-like equation with hamiltonian 
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(fTj) for given quasienergy E and y-component of momen- 
tum k y can be written as a Floquet state: 
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FIG. 2: Transmission probability for central band as a func- 
tion of Qi and 0:3 at <Aq = for two value of 8. 



where the script r — 0, 1,2,3,4 indicates the region, x r 
are the coordinates of the boundaries: xq = x± = —L/2, 
x 2 = —d/2, x$ = d/2, X4 = L/2, (see Fig. 1) and J„ is 
the Bessel function of the first kind. In barrier regions 
(r = l,3): 



s n — s 'n — s gn(E + nhw 

h r - a — J ( E-Vg+nt 
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On the other hand for r = 0, 2, 4: 

s n = s n = sgn(E + nfku), 



h r — h — i ( E + nhu \ _ U2 

K n — K n — y y hvf J K v , 
— 4>n = arctan(fc y /A; rl ) 



(5) 



and J m - n (V r / ftto) — <5 m . n because in this case, the mod- 
ulation amplitude is V r = 0. The phase difference 5 r = 5 
only for r = 3, elsewhere is zero. As Dirac electrons pass 
through a region subjected to time-harmonic potentials, 
transitions from the central band to sidebands (channels) 
at energies E ± mhu (m = 0,1,2,...) occur as electrons 
exchange energy quanta with the oscillating field. 
The wave function is continuous at the boundaries, and 
the continuity condition can be expressed as: 



K n K12 
K21 K22 



(6) 



where the total transfer matrix K = K(0, 1)-K(1, 2) • 
K(2, 3) ■ K(3, 4), and K(r, r + 1) (r = 0, 1, 2, 3) are trans- 
fer matrices that couple the wave function in the r-th 
region to the wave function in the (r + l)-th region. 
We assume an electron propagating from left to right, 
then, A = {<5„,o} and B 4 is the null vector, whereas 
A 4 = {a 4 } and B° = {&^} are the coefficient vectors of 
transmitting waves and of reflecting waves respectively. 
The total transmission probability for quasienergy E is 



T = En=-oo T n, where 



_ S n COs{(j) n ) 4 2 

1 n — ft \ \ a n\ 

S COS{(j)o) 



(7) 



is the probability of the scattering for an electron with 
incident quasienergy E in region into the sideband 
with quasienergy E + nhio in region 4. The mini- 
mum number N of sidebands that need to be consid- 
ered is determined by the strength of the oscillation, 
N > max(Vi/Hui,Vs/Hui) [1, [HJ], and the infinite series 
for T can be truncated to consider a finite number of 
terms starting from —N up to N. 

Numerical calculations have been made for barriers and 
well of the same width: d = 50 nm, L = 150 nm, 
V = 200 meV, to = 5 x 10 12 Hz and A = 50 nm for 
the wavelength of incident electrons. 
Normal incidence. 

The dependence of the transmission probability for the 
central band (To) on ot\ — Vi/fko and = Vs/fko for 
normally incident electrons is shown in Fig. 2 for two 
values of the phase difference S. For small values of a% 
and as the transmission is perfect and take place only 
through the elastic band, because the oscillating barrier 
can be treated approximated by a static one, and for 
normal incidence, perfect transmission (Klein tunneling) 
through static single [l3| and double barriers has been 
obtained [T6l | . With increasing a\ and 0:3, higher and 
lower sidebands become important since electrons can ex- 
change a large number of photons with the time-periodic 
field, decreasing the transmission probability through the 
central band (To) because these probabilities are now 
spread. However, numerical calculations show the ex- 
istence of a critical value 8 C — 2.64 of the phase dif- 
ference, at which the transmission probability is exclu- 
sively through the central band even for large a\ and 03 
whenever \a\ — < 1. Thus, the quasienergy of the 
transmitted electrons is sensitive to the phase difference, 
because tuning ot\ , 013 and 5 could eliminate inelastic side- 
bands doing the function of an energy filter. 
In order to understand the phenomenon behind this re- 
sult, an analytical expression for To at normal incidence is 
derived. Let us take the barriers with the same modula- 
tion amplitudes ai = 0:3. To simplify our calculation and 
get some direct insight into the interference effect [T7| . 
we take ai as a perturbation, and include the lowest or- 
der corrections up to (Vi/hu>) 2 , retaining only the terms 
corresponding to the central band and first sidebands. 
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Hence, the probability for tunneling through the central 
band can be written as: 



To 



7+|(w) : 



(8) 



The term 7 = (zp* + z*p) /8 contains all the interference 
effects. The complex values z and p are given by: 



z — 2gj( l5 +' £ oi2+£i(<?-i+?o+9i)) \^ e iq-iLi _|_ gigiiij 
_|_ e «(fe-ii2+2-Li(g-i+go+gi)) ^ e -2iq^ 1 L 1 _|_ g-2i<j Lij 

_2gi(fc-ii2+£i(9-i+gn+2gi)) 
_|_ e i(2i5+feiL2+2Li(g_i+(;o+gi)) fg— 2iqiLi e ~2iq L 1 
_2 e 4 ( 25 + fe i- L 2+£i(2g-i+go+gi)) 



(9) 



and 



p = e 



i{S+k L 2 +2L 1 {q^ 1 +q 1 )) 



(10) 



Thus, when 7 = 1 the suppression of inelastic side- 
bands is obtained under the condition: 



S c ± (k ±1 - k ) L 2 ± (g - q±i) Li 



(11) 



The left-hand side of (fTTj) is the phase difference between 
two amplitudes. The first one corresponds to particles 
that propagate through the DBS and absorb (emit) the 
energy quantum haj at the first barrier, whereas the sec- 
ond one corresponds to the same propagation with ab- 
sorption (emission) at the second barrier. Noteworthy 
that ± (k±i — ko) L 2 , is the phase difference between the 
amplitude corresponding to particles that absorb (+) or 
emit (— ) a photon close to the left barrier and traverse 
the well with energy E ± hu and the amplitude corre- 
sponding to particles that absorb (emit) a photon close 
to the right barrier and traverse the well with energy E. 
In the same way, ± (go — <7±i) L\ is the spatial phase dif- 
ference between the above mentioned amplitudes through 
the barrier region. Therefore, the suppression of trans- 
mission through the inelastic sidebands is due to destruc- 
tive quantum interference between these two amplitudes. 
At normal incidence k n — fe n _i = q n -i — <ln = w/vf 



and fc_ T 



-(n-1) — <7-(n-l) 



-Lu/vf with 



(n = 1 , 2, ...), result from the gapless and linear electron 
dispersion law near the Fermi energy. Therefore, we can 
also write the expression (If f I) in general form as: 



S c + L 2 + Li = 7r, 



(12) 



where L\ and L 2 are expressed in units of Vf/uj and the 
critical phase difference depends only on barrier and well 
width. Thus, Eq. (fT"2|) indicates a simultaneous destruc- 
tive interference of contiguos channels at <fio — 0. Fig. 
3a) shows the transmission probability through the elas- 
tic channel as a function of L\ and S at ot\ = 7 cal- 
culated after Eq.([7]). The straight line in this figure is 
Ea. (TT2l at L 2 = 0.25. In panel b) only the region where 
I — To < 0.005 for different well widths is shown. The an- 
alytical result presented as condition (|f 2p remains valid 



even in the regime u\ 3> I, when a large number of 
channels coexist, and show excellent agreement with nu- 
merical solutions in Fig. 3. Thus, quantum interference 
between first order process is the physical basis behind 
Eq. (fT2]) . Now, it is clearly seen why \ot\ — a.§\ < 1 is nec- 
essary for the suppression of the transmission through 
inelastic sidebands. The same amplitude of probability 
for the process of absorption (emission) on both barrier 
region is required to achieve destructive interference. It 
should be pointed out that, at normal incidence, perfect 
total transmission, reported for static [IH and oscillat- 
ing [I5| single barriers persists for harmonically driven 
DBS regardless of interference effects. 
Oblique incidence. 

Now, resonant tunneling through DBS at oblique inci- 
dence under time periodic field is investigated. Due to 
the existence of evanescent modes inside the static bar- 
rier [bij ]. the transmission, as a function of the incident 
energy, has a gap which can be tuned changing the height 
of the barrier and/or the angle of incidence [18( . Because 
of coupling of evanescent states in the barrier with res- 
onant states in the well, several resonance peaks with a 
unity value appear in the transmission gap [191 ] . In the 
following, the influence of oscillation amplitudes and in- 
terference effects on resonant tunneling is studied. Fig. 4 
shows total transmission probability as a function of the 
energy of incident electrons at </>o = vr/18 for different 
amplitudes of the oscillating field and 5 = 0. In Fig. 4a) 
the transmission probability for the static double barrier 
is shown as it corresponds to ot\ = 0:3 = 0, which was pre- 
viously obtained. Resonant peaks are narrow and could 
have very important applications in high-speed devices 
based on graphene as has been suggested in Ref. (l9j . 
Considering the effect of PAT, some satellite peaks ap- 
pear on both sides of the two main resonant peaks, and 
as the amplitude of the oscillating field increases, satel- 
lite peaks move away from the two main resonant peaks 
due to emission or absorption of a greater number of pho- 
tons. Coupling between evanescent wave inside the barri- 
ers and propagating wave inside the well occur for several 
energy values. The positions of the satellite peaks are 
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FIG. 3: (Color online) To as a function of L\ and S at (j>o — 0, 
ai = 7 and a) L 2 = 0.25 b) well widths; O.f, 0.25 and 0.4. 
The red line is (JT^J at L 2 fixed. 
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FIG. 4: Total transmission probability on incident energy at 
4>o = 7r/18 and S — for several value of the modulation 
amplitude. 




FIG. 5: To/T as a function of incident energy and angle for 
cti — a 3 = 7 and 8 = 2.64. 



Eh.,2 ± nhjj where are the positions of the main res- 
onant peaks. At a± = ay, = 3 only two main peaks with 
few small satellites are obtained. However, at high am- 
plitude modulation [a% = a% = 7), the two main peaks 
lose its dominance and the intensities of main peaks and 
satellite peaks are reversed. Therefore, strong quenching 
of resonant transmission with increasing amplitude of os- 
cillating field is found as shown in panel d) of Fig. 4. 
For non-normal incidence k± n — k±r n -i) and q±( n -i) ~ 
q± n depend on E ± nhui, Vq and <j>Q. Then, it is not 
possible to determine analytically a critical phase differ- 
ence which cancels simultaneously all contiguos channels. 
However, suppression of inelastic sidebands at S = 2.64 
is obtained numerically not only for normal incidence 
but (To > 0.95 T) for different values of the incident en- 
ergy and incident angle at a.\ = = 7 as shown in 
Fig. 5. Thus, destructive interference effects are robust, 
since the phase diagram (E vs (fro) m Fig- 5 exhibits 




FIG. 6: To /T as a function of phase difference at ai — as — 7. 
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FIG. 7: (Color online) T as a function of incident energy at 
cj>o — n/18 for S — (black line) and 5 — 2.76 (dashed red 
line). 



many white zones at oblique incidence, corresponding 
to strong quenching of the transmission through inelas- 
tic sidebands. These zones change with phase difference 
and become highly anisotropic due to the chiral nature of 
quasiparticles in graphene. Fig. 6 shows Tq/T as a func- 
tion of S for incident energy corresponding to first main 
resonant peaks (E\ — 193.695 meV) and <fio = tt/18. 
Note that transmission is practically only through central 
band for S — 2.76 and for other phase differences inelastic 
sidebands will become more important. Thus, changing 
S we shall observe radical influence in the resonant fea- 
tures. Fig. 7 shows total transmission probability as a 
function of energy at <po = 7r/18 and a\ — = 7 for two 
values of the phase difference. Adjusting the value of 
S resonant transmission can be controlled. For S = 2.76 
satellite peaks are strongly suppressed and almost perfect 
transmission is obtained for E\ = 193.695 meV whereas 
at 5 = a large number of sidebands coexist and res- 
onant tunneling practically disappears. Thus, quantum 
interference makes a considerable contribution to total 
transmission at non-normal incidence. 
Oscillating well. 

In addition, total transmission probability has been also 
calculated when the harmonic potential is applied at the 
well instead of the barrier. In Fig. 8 PAT in a driven DBS 
with oscillating quantum well (panel b) and oscillating 
quantum barriers (panel a), are compared. We observe 
that the driving field with amplitude 3hu> is more effec- 
tive in quenching of resonant tunneling when applied to 
the well. It occurs because, for the incident energy in the 
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FIG. 8: T as a function of incident energy at 0o = 7r/18 for a) 
oscillating quantum barriers and b) oscillating quantum well. 

gap region (Vo — hk y Vf < E < Vo + hk y vj), the evanes- 
cent modes appear inside the barriers and the electrons 
reflect back and forth several rounds in the well. Thus, 
for oscillating quantum well the electrons have enough 
time to interact with the driving field contrary to oscil- 
lating quantum barriers. Consequently, resonant peaks 
are more strongly quenched in driven well. 
In summary, we have carried out a study of electron 
PAT through graphene-based symmetric double barriers 
driven by a periodic potential. Barriers oscillate with the 
same frequency w, different amplitudes and phase differ- 
ence (5. The time-periodic electromagnetic field generates 
additional sidebands at energies E±nhw in the transmis- 



sion probability due to photon absorption or emission. At 
normal incidence, perfect total transmission probability 
(Klein tunneling) persist for harmonically driven DBS. A 
critical phase difference is found such that, total suppres- 
sion of inelastic sidebands due to destructive interference 
between waves is obtained. A condition for the simul- 
taneous cancellation of contiguos channels as a result of 
linear dispersion law, is derived for small amplitude mod- 
ulation and is valid even in the regime (ai = Q13 > 1). 
Thus, energy of the transmitted electrons is sensitive to 
phase difference and inelastic sidebands can be removed 
by tuning structural parameters of the DBS according 
to Eq. (fl"2|) . For oblique incidence, suppression of inelas- 
tic sidebands occurs and destructive interference plays 
a fundamental role in the total transmission probabil- 
ity. Resonant transmission may be regulated, as resonant 
tunneling is quenched at 5 — and recovered for cer- 
tain phase difference. Moreover, quenching of resonant 
transmission in driven well is more drastic than driven 
barriers. In conclusion, we have shown that quantum in- 
terference has an important effect on quasiparticles tun- 
neling through a time-dependent graphene-based double 
barrier. This phenomenon has potential applications in 
graphene-based electronic devices such as energy filters 
and high-frequency radiation detectors. 
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